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We apply CMB lensing techniques to large scale structure and solve for the 3-D cosmic tidal field. 
We use small scale filamentary structures to solve for the large scale tidal shear and gravitational 
potential. By comparing this to the redshift space density field, one can measure the gravitational 
growth factor on large scales without cosmic variance. This potentially enables accurate measure- 
ments of neutrino masses and reconstruction of radial modes lost in 21 cm intensity mapping, which 
are essential for CMB and other cross correlations. We relate the tidal fields to the squeezed limit 
bispectrum, and present initial results from simulations and data from the SDSS. 



Introduction. - The large scale structure of the uni- 
verse shows striking non-Gaussian features, often de- 
scribed as the cosmic web. The filamentary nature arises 
from gravitational tidal shear. It is the same principle 
that leads to ocean tides on earth: the residual anisotropy 
of gravitational forces in a free falling frame. The strong 
non-Gaussian nature of this system has traditionally led 
to a reduction in cosmological information that can be ex- 
tracted from a large survey [HE]. In this paper, we turn 
the process around, and exploit the tidal non-Gaussianity 
to improve the measurement of large scale structures. 

The gravitationally induced displacement is a three 
component vector field. Only changes in displacement 
are observable, which are described by the Jacobian of 
the displacement field. The Jacobian can result in a 
change of volume and of shape. In this paper we shall 
study the change of shape, which is more robust since 
many other processes can lead to a change of number 
density. This approach is equivalent to the gravity wave 
shear computed in [3] . We are using small scale structure 
alignment to solve for the large scale tidal field. This pro- 
vides many independent samples to accurately measure 
the large scale tidal field. 

Generally speaking, the shape tensor in three dimen- 
sions is described by 5 numbers: three Euler angles and 
two axis ratios. The gravitational potential is a single 
number, which is five fold overdetermined. The radial 
changes are affected by peculiar velocities, which are be- 
yond the scope of this paper, hence we only use the shear 
in the plane of the sky. This tangential shear is described 
by two numbers, and we will use the gravitational lens- 
ing notation, where the two variables are called 71, 72 [IJ. 
A linear transformation decomposes them into a diver- 
gence, or E-like mode (called k), and a curl, or B-like 
mode[SHZ]- 

Tidal Reconstruction. - The reconstruction of tidal 
shear is described by the same formulation as the re- 
construction of gravitational lensing induced shear. It 
leads to a local anisotropy of quadratic statistics. As 



shown in jS], there exists an optimal quadratic estima- 
tor that is solvable but distinct for both Gaussian and 
Non-Gaussian fields. In the first case, the optimal esti- 
mator can be expressed in terms of the power spectrum of 
the density field alone. In general, however, the optimal 
weights need to be computed from simulations [3]. 

Heuristic kernel. - Here we work through a simple, 
slightly sub-optimal, scenario. We start with the density 
field 5(x). In analogy to CMB lensing, quadratic esti- 
mators are outer products of gradients. The first step is 
a convolution, which filters for the small scale structure 
(the gradient downweights the large scales). The large 
scale gravitational field <f> will be a linear convolution 
of this local small scale quadratic estimator. For sim- 
plicity, we use a Gaussian window W(r) = exp J. 

We obtain a smoothed density field: 8{x) = J W{x — 
x'^x'^x'. 

Because the quadratic estimator heavily weights the 
high density regions, we further apply a Gaussianization 
technique [TU] to the non-Gaussian (smoothed) density 
field by taking 8 g = log(l + 5). Following gravitational 
lensing procedures, we construct the two shear compo- 
nents as 71 = (d x 6 g ) 2 - (dySg) 2 and 72 = 2(d x S g )d y S g . 
We can then reconstruct the dark matter field k, in anal- 
ogy with the convergence field, which is a linear convo- 
lution of these two estimators. In terms of differential 
operators, we define d = (d x — dy)ji + 2d x d y j2, and 
solve Poisson's equation for the density (d x + dy)K z = Nd 
on each z slice. The 3-D dark matter density field 
n(x,y, z) is given by one more ratio of wave numbers 
(d x + dy)n = (d x + dy + d z )n z . N is a normalization 
constant which is derived in Lu and Pen [2]. 

Simulations. - We ran N-body simulations with the 
CUBEP3M code[TT], evolving 256 3 particles on a 512 3 
grid in a 322/i _1 Mpc box. The smoothing window W 
has a width a = 1.25ft. _1 Mpc. Figure [l] shows a slice of 
the original and reconstructed smoothed density fields, 
both once again smoothed on a 8/i _1 window to reduce 
the small scale noise. We found little dependence on the 




FIG. 1: Density field slice smoothed on 8/i _1 Mpc in a 
322/i -1 Mpc simulation box. Top panel: Original dark matter 
field. Bottom panel: reconstructed smoothed field. 



windowing scale until we smooth on linear scales, which 
is related to the information saturation phenomenon dis- 
cussed below. 

Figure [2] shows the raw and reconstructed power spec- 
tra, as well as the cross correlation spectrum. The ampli- 
tude of the reconstructed spectrum was scaled to match 
that of the dark matter. In principle, this amplitude is 
determined for truly Gaussian random fields, but in prac- 
tice these assumptions do not hold precisely. The am- 
plitude of the cross correlation is now completely deter- 
mined, and we see a good cross correlation over a decade 
in wave number. 
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FIG. 2: Simulation results: Dark matter power spectrum 
(solid line), tidal reconstructed power spectrum (dashed line), 
cross power spectrum (dotted line). We see good agreement 
between the reconstructed power spectrum and the dark mat- 
ter power. 

Data. - We also applied the technique to a volume lim- 
ited subsample of the Sloan Digital Sky Survey (SDSS) 
in which 112043 galaxies were selected. We constructed 
a Cartesian mapping as follows: A 0.9ft- -1 Gpc box is 
centered on the north galactic pole, at z = 0.1. We used 
the redshift range 0.7 < z < 1.3. The angular position is 
mapped at SOOh" 1 Mpc per radian, and the radial posi- 
tion at O.Ol/i -1 Mpc/ (km/sec). A random catalog with 
five times as many galaxies was used as a density nor- 
malization. The log of the smoothed random density is 
subtracted from the log of the smoothed galaxy count. 
Pixels with less than half the peak random galaxy den- 
sity count were masked. Figure [3] shows a slice of the 
galaxy and reconstructed density fields at z = 0.1. 

The original and reconstructed power spectra are 
shown in Figure [4j where we have again scaled the re- 
constructed power spectrum. When comparing the bi- 
ased galaxies, the bias will always be a free parameter, 
such that the non-linear effects of normalization can be 
absorbed as a relative bias. The cross power spectrum 
has a cross correlation coefficient r > 0.5 for all k up to 
0.2/i/Mpc, thereby demomstrating a qualitative success- 
ful dark matter reconstruction. 

Discussion. - We have seen a successful application of 
a heuristic tidal field estimator. Here we will discuss the 
theoretical framework, and future directions. 

The large scale tidal shear field is the expectation value 
of the quadratic outer product 

Tij ex ((cV)cy). (1) 




FIG. 3: z = 0.1 slices of SDSS survey embedded in a 900ft" 1 
Mpc box, smoothed by a window function. Top panel: Orig- 
inal galaxy field. Bottom panel: dark matter map recon- 
structed from small scale galaxy density field alignments. 



Redshift space distortions complicate every components 
that include i = 3. Our shear components are related as 
7i = Tn —T22 and 72 = 2Ti2. It was shown that this two 
index tensor is proportional to the traceless tidal shear 
tensor (see eqn(35) of [2]) derived from the gravitational 
potential: Ty = didj<j) — ^-V 2 0. The trace, which corre- 
sponds to the local mean density, in principle transforms 
in a similar way, however we will defer such study to the 
future. The trace free estimator is quadratic in the den- 
sity field, and does not require second order perturbation 
theory to compute. Being trace free, it does not backre- 



FIG. 4: SDSS results: galaxy power spectrum (solid line), 
tidal reconstructed power spectrum (dashed line), cross power 
spectrum (dotted line). We see good agreement for the recon- 
structed power spectrum. 

act on the mean density. We are only considering waves 
in <j) on much larger scales k^ than the galaxy density 8: 
k s > fc^. 

The reconstructed density field is derived from a con- 
volution: 

k = J Ti^x^Kijix - x')d 3 x' (2) 
where the kernel is given in Fourier space as 

In the limit that the tidal mode is linear and has a 
wave length much longer then the smoothing scale of the 
Gaussian window, the perturbative coupling arises solely 
from the induced gravitational displacement. 

Free falling observers are only sensitive to tidal forces, 
and the tidal force on the window scale is just a stretching 
of coordinates. The tidal distortion field can be viewed as 
an integral over the history of gravity. We note that the 
actual density field only follows linear evolution weakly, 
even on apparently linear scales [12]. 

Relation to CMB Lensing. - The effect described here 
is identical to the CMB lensing phenomenon, which also 
induces a three point correlation between two points on 
the CMB and one in the intervening space [13]. In our 
case, two wave vectors are on small scales, and one on 
large scales, all in the same volume. 

Connection to Non-Gaussianity. - Our separation of 
scales is equivalent to bispectrum calculations in the 
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squeezed limit. We are measuring the variance modu- 
lated by a large scale, small k mode, again in analogy 
with the CMB case [14"] . The coordinate change is the 
leading order effect relevant for the shear, whose mea- 
surement corresponds to using only the term dependent 
on the angle a between the small k$ vector and the two 
large ones ks- Figure [4] is the covariance of (5), and k 
is in turn a quadratic function of 8 at large wave num- 
bers. We have effectively integrated the bispectrum over 
its quadrupole f (5(k r f > )\S(kg)\ 2 ) cos(2a)da. 

This is the leading effect contributing to the informa- 
tion saturation phenomenon [T]. It was found that the 
variance in the power spectrum on quasi-linear scales was 
up to three orders of magnitude larger than expected for 
Gaussian random fields. This was also observered for 
variances in the shear estimators [51 [S]. It had been pro- 
posed that this saturation was due to Poisson noise in 
the virialized halos[15], however Poisson noise does not 
lead to an increase in the shear variance, hence can not 
be the complete picture. The shear has been numerically 
observed to saturate at a similar information content as 
the mean variance. Furthermore, reconstruction tech- 
niques are able to increase the information content [IB] . 
which is not expected in a Poisson model. 

The cosmic tides picture can qualitatively explain this 
information saturation effect. A 10% tidal distortion re- 
sults in a systematic change of the variance by ~ 10%. 
When more than 100 modes are measured, their vari- 
ance is dominated by the larger scale shear. An analo- 
gous effect applies to the power spectrum variances. This 
implies that saturation is dominated by modes near the 
peak of the power spectrum, k ~ 0.02. 

Recovering Lost 21cm Modes. - 21cm intensity map- 
ping has emerged as a promising technique to map the 
large scale structure of the universe, at redshifts z from 1 
to 10. Unfortunately, many of the key cross correlations 
with photo-x galaxy and the CMB have been thought to 
be impossible due to foreground contamination for radial 
modes with small wave numbers[T71 [TH]. Our proposed 
tidal reconstruction technique works best in this regime, 
and opens up a new set of possibilities. 

Reducing Sample Variance. — The high cross correla- 
tion seen in the simulation suggests that one can measure 
the dark matter density on large scales without redshift 
space distortions. By comparing the galaxy power spec- 
trum to the k power spectrum as a function of angle to 
the line of sight, one can solve for the velocity contribu- 
tion. This allows for a measurement of velocity on the 
same mode as the density, hence determines the veloc- 
ity growth factor without sample variance, analogous to 
McDonald and Seljak [19]. In principle, this could enable 
precision measurements of neutrino masses. 

Conclusion. - We have presented a new framework to 
measure the effects of gravity through its tidal effect on 
the structure of galaxy clustering. We isolated this effect 
theoretically from other complications by projecting out 



mean density changes. The effect becomes analogous to 
the lensing shear on large scale structure or the CMB. 
We have tested the dark matter density reconstruction 
in simulations and with the SDSS data. We found good 
agreement, with reconstruction noise less than sample 
variance. This opens up the window for precision mea- 
surements of the transfer function, potentially measuring 
neutrino masses or testing modified gravity. 
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